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Abstract 

We argue that due to various restrictions cosmic strings and monopole- 
string networks are not likely to produce the observed flux of ultra-high en- 
ergy cosmic rays (UHECR). Among the topological defects studied so far, the 
most promising UHECR sources are necklaces and monopolonia. Other viable 
sources which are similar to topological defects are relic superheavy particles. 
All these sources have an excess of pions (and thus photons) over nucleons at 
production. We demonstrate that in the case of necklaces the diffuse proton 
flux can be larger than photon flux, due to absorption of the latter on ra- 
diobackground, while monopolonia and relic particles are concentrated in the 
Galactic halo, and the photon flux dominates. Another signature of the latter 
sources is anisotropy imposed by asymmetric position of the sun in the Galac- 
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tic halo. In all cases considered so far, including necklaces, photons must be 
present in ultra-high energy radiation observed from topological defects, and 
experimental discrimination between photon-induced and proton-induced ex- 
tensive air showers can give a clue to the origin of ultra-high energy cosmic 
rays. 

I. INTRODUCTION 

The observation of cosmic ray particles with energies higher than 10 11 GeV |3J gives 
a serious challenge to the known mechanisms of acceleration. The shock acceleration in 
various astrophysical objects typically gives maximal energy of accelerated protons less than 
(1 — 3) • 10 10 GeV @ (see however ||). The unipolar induction can provide the maximum 
energy 1 ■ 10 11 GeV only for extreme values of the parameters Q]. Much attention has 
recently been given to acceleration by ultrarelativistic shocks [[J, ||. The particles here 
can gain a tremendous increase in energy, equal to T 2 , at a single reflection, where F is 
the Lorentz factor of the shock. However, it is known (see e.g. the simulation for pulsar 
relativistic wind in f7|) that particles entering the shock region are captured there or at least 
have a small probability to escape. 

Topological defects (TD) (for a review see ||) can naturally produce particles of ultrahigh 
energies (UHE) well in excess of those observed in cosmic rays (CR). In most cases the 
problem with topological defects is not the maximum energy but the fluxes. However, in 
some cases the predicted fluxes are comparable with observations. 

Usually, UHE particles appear at the decays of superheavy (SH) particles produced 
by TD. (We shall refer to these SH particles as X-particles). Examples discussed in 
the literature include ejection of X-particles from superconducting strings, emission of X- 
particles from cusps or intersections of "ordinary" strings, and production of such particles 
in monopole-antimonopole annihilations. Metastable SH particles can also be relics of an 
earlier epoch, produced by a thermal or some other mechanism in the early Universe. 
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A rather exceptional mechanism of UHE particle production is given by radiation of 
accelerated monopoles connected by strings. In this case a monopole can radiate gluons with 
very large Lorentz factors and with virtualities of the order of the monopole acceleration. 

A common signature of all extragalactic UHECR is the Greisen-Zatsepin-Kuzmin (GZK) 
cutoff |J. It reveals itself as a steepening of the spectrum of UHE protons and nuclei due 
to their interaction with microwave radiation. The steepening starts at E « 3 • 10 10 GeV. 
Apart from this steepening there is another signature of interaction of extragalactic CR 
with microwave radiation: a bump in the spectrum preceeding the cutoff. The bump is a 
consequence of the proton number conservation in the spectrum: protons loose energy and 
are accumulated before the cutoff. 

In this paper we will discuss the signatures of UHECR from TD distinguishing them 
from particles produced by astrophysical accelerators. 

We will confine ourselves here to the case of the conventional primary particles, protons 
and photons, and will not consider the other UHE signal carriers discussed in the literature 
such as neutrinos pIO], Lightest Supersymmetric Particles JT^JTBI], relativistic monopoles 



TA \ and vortons [I5|. 



Throughout the paper we shall use the following numerical values and abbreviations: the 
dimensionless Hubble constant h = 0.65, the Cold Dark Matter density in terms of critical 
density Qcdm = 0-"2h 2 , the size of Dark Matter halo R h = 100 kpc, UHECR - for Ultra 
High Energy Cosmic Rays and UHE - for Ultra High Energy, TD - for Topological Defect, 
SH - for Superheavy, CDM - for Cold Dark Matter, SUSY - for Supersymmetry, LLA - for 
Leading Logarithmic Approximation, AGN - for Active Galactic Nucleus, GC and AC - for 
Galactic Center and Anticenter, respectively. 

II. CONSTRAINTS AND SIGNATURES 

A common characteristic feature of UHE particle production in TD is an excess of pions 
over nucleons PIHTf • As a result in many cases the observed UHE gamma-ray flux dominates 
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over proton flux and it makes 7/p ratio "a diagnostic tool" [16] for TD. We shall discuss 
this signature quantitatively for various sources, such as cosmic necklaces, monopolonia and 
superheavy relic particles. 

An excess of photons over protons at generation is present in all known mechanisms 
of high energy particle generation: decay of SH particles produced by strings and cusps, 
annihilation of monopoles and radiation of UHE particles by accelerated monopoles. The 
order of magnitude of this excess can be estimated from the ratio of energy transferred to 
pions and nucleons in the QCD cascade. For example, N/ n ratio for the decay Z° — > hadrons 
is about 5% (N includes p, p, n, n and ir - charged and neutral pions). To estimate 7/A-ratio 
at energy of interest, one should take into account two photons produced by each 7r°-decay, 
the fraction (1/3) of neutral pions relative to all pions and the energy spectrum of produced 
hadrons. It gives roughly 7/iV ~ 10 at production. However, in the observed diffuse flux 
at E > 10 6 GeV the proton component [] can dominate because of the strong absorption 
of high energy photons on the background photon radiation. UHE photons are absorbed 
mainly on radio background [|TR|ni] and the absorption length is sensitive to low-frequency 



cutoff in this background. 

The pion dominance of hadron production by TD has two consequences, (i) A large part 
of hadron energy is transferred, due to pion decays, to an electromagnetic cascade. This 
can be used to derive an upper bound on the UHE proton flux, (ii) TD localized inside the 
sphere of gamma-ray absorption give a direct flux of UHE photons at the Earth, producing 
thus an observable anisotropy. 

The e-m cascade upper limit on UHECR arises due to cascading of electrons and photons 
in the Universe down to the observed energies 10 MeV — 100 GeV. The flux of the cascade 
photons at these energies must be lower than the extragalactic flux measured by EGRET 
pOfl. The cascade photon flux is below the EGRET extragalactic flux at 10 MeV— 100 GeV if 



1 Here and below we shall refer to the p + p diffuse flux from TD as to proton flux. 
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the energy density of the cascade photons is u cas < 2 • eV/cm 3 . This result follows from 
Monte Carlo simulation in pl| , p2| and from our own estimates based on analytic calculations 
in Pp3|). This limit (which will be used in our calculations) is rather rigorous: due to 
uncertainties in infra-red flux and intergalactic magnetic field the allowed density of the 
cascade radiation can be as high as (3 — 5) • 10~ 6 eV/cm 3 . 

The energy density of e-m cascades can be readily calculated as 

1 r f o 



UJcas = -^firmx \ dtn x (t) (1 + z) , (1) 







where hx{t) is the rate of X-particle production at the epoch t (redshift z), to = 2.06 • 
10 17 /i _1 s is the age of the Universe, h is the dimensionless Hubble constant, f n is the 
fraction of energy transferred to pions at the decay of X-particle, and 1/2 takes into account 
that half of this energy goes into e-m cascade. One can parametrize the effect of the evolution 



of TD on X-particle production as nx(t) = nx(to)(t/to)~ m . In most cases, e.g. ordinary 
strings, monopolonium, and necklaces, m = 3. In this case the integral in Eq.(|I]) reduces to 
/ dz/(l + z) 2 , i.e. the evolutionary effects are absent. In the case of superconducting strings, 
m is model-dependent and it can be m = 4 or larger. Weak cosmological effects are present 
in this case. 

With u cas from Eq. (1) and the diffuse proton flux being determined by /atWx(^o), where 
/tv is a fraction of mx transferred to UHE nucleons, one can obtain an upper bound on the 



diffuse flux of UHE (for another approach see ) . 

It is easy to generalize the calculations above to the case when X-particles are produced 
with a large Lorentz factor V . 

Another general restriction on TD models is given by the distance between the topological 
defects, D. There are three distance scales in our problem: the distance between TD, 
D, the photon absorption length, R^(E), and the proton attenuation length, R p (E) = 
ciE^dE/dt)- 1 , where dE/dt is the energy loss of UHE proton on microwave radiation. We 
shall analyze the case _R 7 < R p , though at very high energies, E > 10 12 GeV, they can be 
comparable. 



For TD with D > R P (E) the diffuse flux at a representative point between TD is ex- 
ponentially suppressed. For a power-law generation spectrum with exponent 7 g and the 
distance D/2 to the nearest source, the suppression factor for rectilinear propagation is 

exp (-( % -l)D/2R p (E)). (2) 

In the case of diffusive propagation, D/2 should be replaced by the propagation time. Such 
TD are disfavored as sources of UHECR, because their spectrum either has an exponential 
cutoff at energy where R P {E) < D or, in case of accidental proximity of a source, the flux 
is anisotropic. The anisotropy can be estimated as 5 ~ l/2ct, where I is the distance to the 
source and t is the propagation time. 

In the other extreme case, D < R 1 (E) 1 a TD located inside the photonic sphere R 1 
creates a direct UHE photon flux, F 7 = Q^/inr 2 at the point of observation. The produced 
Extensive Atmospheric Showers (EAS) can be identified as photon-induced EAS in the 
direction of the source. 

The proton-induced showers can dominate in the case R^(E) < D < R P {E). The sources 
might be absent inside the photonic sphere with radius i? 7 and thus, if magnetic field is 
strong enough, no source is seen directly. The proton spectrum exhibits the usual GZK 
cutoff, due to energy losses of the nucleons produced by the sources beyond the distance R p . 
Since at extremely high energies R p and i? 7 are not much different, the number of sources 
inside the GZK sphere of radius R p is not large and some anisotropy is expected. 

Finally let us turn to the case D « R 1 {E). Generically, this is a case of uniformly 
distributed sources. The proton showers dominate when R P (E) > R^(E)\ their spectrum 
has the usual GZK-cutoff. A certain fraction of showers, ~ R^/Rp, are the photon-induced 
ones, and they correspond to the direct arrival of the photons. The 7/p-ratio depends thus 
on the calculations of attenuation lengths of photons and protons. We shall analyze this 
case quantitatively. 

A special situation arises when the sources are concentrated in the galactic halo. This 
happens in the case of relic superheavy particles |25]. Here we point out that the same 
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phenomenon of galactic enhancement occurs for monopolonium and decaying vortons. 

The energy losses and absorption are negligible for the halo model and thus photons 
strongly dominate. 

Another signature of this model was indicated recently in P6[ : because of asymmetric 



position of the Sun in the Galactic halo, there is a considerable anisotropy of UHE photon 
flux. 

III. FLUXES 

In this section we shall give the formulae for UHE proton and photon fluxes from extra- 
galactic space and from the halo (for the case of SH relic particles and monopolonia). The 
basic quantity which determines these fluxes is the rate of X-particles production fix- In 
the case of an extragalactic flux we assume that sources are uniformly distributed in space 
and rix does not depend on time. For the galactic model we assume that hx is a function 
of the distance from Galactic Center R. 

X-particles can be produced by TD at rest or with a Lorentz factor V. The decay of an X- 
particle results in a parton cascade. The energy spectrum of hadrons outside the confinement 
radius is described by fragmentation functions Wx(x,m x ) for nucleons and ^(^mj) for 
pions, where x = 2E/mx and E is the energy of a proton or a photon. For fragmentation 



functions we use the supersymmetric generalization [£8] of the LLA limiting spectrum of 



QCD cascade |27| normalized by the fraction of energy transferred to the nucleons fx, and 



pions /tt, respectively: 

/ dxxWi(x,m x ) = 2fi (3) 
J o 

with i = N, re. In all calculations below we shall use ~ 0.5, as suggested by calculations 



which show that about half of energy of SUSY-QCD cascade is taken away by neutrali- 
nos and high-energy leptons. For the ratio fx/ 'fir we shall use 0.05 inspired by data on Z°- 
decay. The SUSY-QCD fragmentation function considerably differs from that of ordinary 
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QCD ||28|| : the maximum of the Gaussian peak is shifted towards smaller x, and the peak is 



much higher, while at larger x SUSY-QCD spectrum is below the ordinary QCD spectrum. 
Since u cas is determined by large x, there are more UHE particles at fixed u cas in case of 
SUSY-QCD fragmentation function, as compared with ordinary QCD spectrum. 

As it is discussed in [fZT|-|25| the normalization of analytic solutions by conservation of 
total energy, results in large uncertainty in the normalization constant. On the other hand, 
the ordinary QCD spectrum (used e.g. in f29j) has the the shape of the spectrum much 
different from SUSY-QCD spectrum. One needs SUSY-QCD Monte Carlo simulation for 
properly normalized SUSY-QCD spectrum. 

The spectra of extragalactic UHE protons and photons can be calculated as 

I p (E)-—J o dtW N (x g ,m x )^^, (4) 



lf r (E) = Ry(E)—^- / ^WA^mx), (5) 

imi X J2E/m x x 



dE g (E,z g ) _ n ( ( [* 3 dz_ 



(]E (1 + z 9 )exp [J 3 ^ + z) 1/2 (db(E, 0)/dE) E=Eg{z) j , (6) 

where x g = 2E g /mx and E g (E,t) is the energy of generation at time t for a proton with 
energy E now. The energy E g is determined by proton energy losses on microwave radiation 
dE/dt = b(E, z) and due to redshift. The expression dE g /dE, given by Eq.(||), is taken from 

The photon absorption length, i? 7 (i?), at very high energies is determined mostly by pair 
production on radio background. For our calculations we use the absorption lengths from 
]Tj| . The photon absorption lengths from [H| and [19] are plotted in Fig.l. 



As one can see from Eq.(|4]), the proton flux is calculated with recoil protons taken 
into account, while for photons we neglect multiplication due to cascade on background 
photons. The reason is that electron-positron pairs produced as a result of absorption of 
primary photons on radio-photons, are losing their energies on radiobackground practically 
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continuously and most of their energy is lost before a collision with a microwave photon - 
the process responsible for the cascade development. 

The ratios j/p = I 7 (E)/I p (E) as determined by Eq's((|, |5]) are plotted in Fig.2 as function 
of energy for two extreme cases of absorption length from . One can see that even for 
the exceptional case of low j/p ratio considered here (D <C i? 7 ), it becomes appreciable, 
0.5 < 7/p < 2.5, at energy 3 • 10 11 GeV and increases further with energy. 

Let us consider now the halo model. In cases of SH relic particles |25|] and monopolonia, 



their density is everywhere proportional to the density of Cold Dark Matter (CDM). Thus, 
the ratio of the rates of X-particle production in the halo, n x (r Q ), and in extragalactic 
space, h e x tr , is given by 

r^extr O^.' ^ ' 

n X i t'CDMPcr 

where r Q is the distance between the Sun and the Galactic Center, VL C dm is the CDM density 
in the extragalactic space in units of critical density p cr and Pcdm^J'®) ~ 0-3 GeV/ cm 3 is 
the local CDM density in the halo. 

The fluxes of UHE protons (i = p) and photons (i = 7) from the halo can be calculated 



as 



iHE) = / cos Od cos 6 / drh h x (R), (8) 

4:71 JO JO 

where 9 is an angle relative to the direction of GC and 



r m ax (0) = r Q cos 9 + J R\ - r% sin 2 9. 



The rate of particle production, h x (R), is parametrize here as density of DM [31 



Ux ^ (R/r p[l + (R/r ) a ]il 3 -ry^ (9) 

where Hq is the normalizing rate, R is the distance from the Galactic Center, a, (3, 7 = (2, 2, 0) 
correspond to isothermal profile, a, (3, 7 = (2,3,0.2) gives, according to |3T]] , the best fit to 
observational data and a,/?, 7 = (1,3, 1) describes the numerical simulations of Ref.( |32| ). 



The functions W P (E) and W 1 (E) in Eq.(||) can be given in terms of fragmentation func- 
tions, WN^x.mx) and W^o^x.mx), as 

W P (E) = — W N (x,m x ), 
m x 

2 r 1 dr 

W 1 {E) = / —WA^m x ). (10) 

mx J2E/m x X 

Relic particles and monopolonia fill also extragalactic space with space density propor- 
tional to flcDMPcr- The decay rate h e x tr is given by Eq.(^) and the fluxes by Eq's(f|) and 

(!)• 

IV. SOURCES 

The following Topological Defects have been discussed as potential sources of UHE par- 
ticles. 

(i) Superconducting strings [j33fl , 



ii) Ordinary strings [pi], including the cusp radiation [35 



(iii) Networks of monopoles connected by strings [|36|,j3"7 



(iv) Necklaces pSfl : hybrid topological defects, where each monopole is attached to two 
strings, 

(v) Magnetic monopoles, or more precisely bound monopole-antimonopole pairs 
(monopolonium 



(vi) Vortons ||41|| : small loops of superconducting string stabilized by their angular momen- 
tum. 

Finally we include in this list SH quasistable relic particles produced in the early Universe 
n, |25|^3| - ^5l , p9| . Here we shall apply the criteria discussed in the previous section to each 
of these sources. 
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(i) Superconducting strings 



Superconducting strings produce SH particles when the electric current in the strings 
reaches the critical value, i = i c . In some scenarios, e.g. [[46] where the current is induced 
by primordial magnetic field, the critical current produces strong magnetic field, in which 



all high energy particles degrade catastrophically in energy ||47|| . However, for ac currents 
there are portions of the string with large electric charge and small current. High energy 
particles can escape from there. 

Large ac currents can be induced in string loops as they oscillate in galactic or extragalac- 
tic magnetic fields. Even if the string current is typically well below critical, super-critical 
currents can be reached in the vicinity of cusps, where the string shrinks by a large factor 
and density of charge carriers is greatly enhanced. In this case, SH particles are emitted 
with large Lorentz factors. 

Loops can also acquire dc currents at the time of formation, when they are chopped off 
the infinite strings. As the loops lose their energy by gravitational radiation, they shrink, 
the dc currents grow, and eventually become overcritical. There could be a variety of as- 
trophysical mechanisms for excitation of the electric current in superconducting strings, but 
for all mechanisms considered so far the flux of UHE particles is smaller than the observed 



flux 48 . However, the number of possibilities to be explored here is very large, and more 



work is needed to reach a definitive conclusion. 



(ii) Ordinary strings 



There are several mechanisms by which ordinary strings can produce UHE particles. 
For a special choice of initial conditions, an ordinary loop can collapse to a double line, 
releasing its total energy in the form of X-particles However, as noted in ||34|| , the 



probability of this mode of collapse is extremely small, and its contribution to the overall 
flux of UHE particles is negligible. 
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String loops can also produce X-particles when they self- intersect (e.g. p9|). Each in- 
tersection, however, gives only a few particles, and the corresponding flux is very small 



Superheavy particles with large Lorentz factors can be produced in the annihilation of 
cusps, when the two cusp segments overlap [P5fl . The energy released in a single cusp event 
can be quite large, but again, the resulting flux of UHE particles is too small to account for 
the observations |55|^ . 



One effect which was not considered in |50| and []35[ could increase the production of 



UHE particles [f^plfl . As a non- intersecting closed loop oscillates and radiates away its 



energy, the loop configuration is gradually changing. After the loop has lost a substantial 
part of its energy, it is likely to self-intersect and fragment into several smaller loops. These 
daughter loops will go through the same cycle, and the process will continue until the size 
of the fragments becomes comparable to the string thickness, at which point the fragment 
loops disintegrate into relativistic particles. UHE radiation is also emitted from cusps on 
the daughter loops. 

The process of loop fragmentation is not well understood. We do not know, for example, 
the number and size distribution of the fragments. The only numerical simulation of the 



fragmentation process that we are aware of, |p2| , used initial loops of arbitrary shape, and 



it is not at all clear that the results are relevant for loops produced in a realistic evolving 
network. To address this problem, we used a simple analytic model (see Appendix A) which 
assumes that, in each round of fragmentation, a loop looses a fixed fraction (1 — fi) of its 
energy to gravitational waves and breaks into a fixed number (N% + 1) of daughter loops of 
roughly equal size. The daughter loops move with Lorentz factors I\ in the rest frame of 



the parent loop. This model, which is similar to that introduced in Ref. ||51|| , is analyzed 
in Appendix A, with the conclusion that, for reasonable values of the parameters fx, Ni 
and Ti, the UHECR flux from fragmenting loops is still too small. While these results are 
suggestive, a detailed numerical simulation of loop fragmentation with gravitational back- 
reaction will be needed to reach a more definitive conclusion. However, there is one general 

12 



argument against this model of UHECR production. In the chain of fragmentation the 
Lorentz factor of the daughter loops increases multiplicatively, e.g. in our model as V™, and 
at the moment of collapse, reaches a large value. The large maximum energy of X-particle 
results, for the observed UHECR flux, in an unacceptably high cascade gamma radiation. 
This fact, observed in all calculations (see for example pifl ) is easy to understand. The total 
energy of decaying X-particle is released mostly in the particles of highest energies. When 
one fixes the flux of UHECR at a given energy, e.g. E = 1.0 • 10 20 eV, and increases E max , 
the energy transferred to e-m cascades increases. This argument works especially strongly 
in case of cusps, where the Lorentz factor can reach tremendous values, e.g. T ~ 10 15 — 10 17 



in |35|. 



It has been recently argued |k| that long strings lose most of their energy not by produc- 



tion of closed loops, as it is generally believed, but by direct emission of heavy X-particles. 
If correct, this claim will change dramatically the standard picture of string evolution. It 
has been also suggested that the decay products of particles produced in this way can ex- 
plain the observed flux of UHE cosmic rays pB|Jo^ ]. However, we are not convinced that the 
conclusions of Ref. E3[ are justified. In fact, we believe that numerical simulations described 



in allow an alternative interpretation. The initial string separation in these simulations 
is comparable to the string thickness. As a result, string intersections and reconnections 
can generate a large amount of energy in the form of short-wavelength perturbations on the 
strings. These perturbations may then be released in the form of particles from portions of 
the string that undergo contraction in the course of the following evolution. 



Even if the conclusions of |53| were correct, the particle production mechanism suggested 
in that paper cannot explain the observed flux of UHE particles. If particles are emitted 
directly from long strings, then the distance between UHE particle sources is of the order of 
the Hubble distance, D ~ t ^> R p . According to the discussion in the preceeding section, 
the flux of UHE particles in this case is exponentially suppressed, or, in the case of accidental 
proximity of a string to the observer, the flux is strongly anisotropic. A fine-tuning in the 
position of the observer is needed to reconcile both requirements, because long strings are 
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separated by a Hubble distance. 



(in) Network of monopoles connected by strings. 



The sequence of phase transitions 

G —> H x U(l) —> H x Z N (11) 

results in the formation of monopole-string networks in which each monopole is attached to 
N strings. Most of the monopoles and most of the strings belong to one infinite network. 
The evolution of networks is expected to be scale-invariant with a characteristic scale 

d = Kb, (12) 

where k = const. The scale d gives the average distance between monopoles and the typical 
length of string segments. Each string attached to a monopole pulls it with a force equal to 
the string tension, \i ~ rj^, where n s is the symmetry breaking vev of strings. The monopoles 
are accelerated and radiate gauge quanta at the rate 

dE h 2 2 fi 2 



(13) 



dt Qn g 2 m 2 ' 

where h ~ 2ir/g and m are the monopole charge and mass, respectively, g is the gauge 
coupling and a ~ fi/m is the monopole acceleration. The GUT monopole has magnetic 
and chromomagnetic charges, h m ~ 2ii/e and h s ~ 2ir/3g s , respectively, where e is the e-m 
coupling constant and g s is the color coupling constant. From Eq.flTB]) it follows that the 
energy losses are dominated by e-m radiation. Then a simple energy balance analysis gives 
the value of k in Eq.(|i~2"|) f36|j , 



k ~ fi/e m 2 . (14) 

The energy of gauge quanta (practically photons and gluons) radiated by monopoles can be 
estimated assuming a rough energy equipartition between the monopole and string subsys- 
tems. Then the monopoles have a typical energy E ~ fid and Lorentz factor T]\j ~ fid/m. 
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If the mass of gauge quanta (or the virtuality Q 2 in the case of gluon) is smaller than the 
monopole acceleration a, the typical energy of gauge quanta is e ~ T M a [|57|| ; otherwise the 
production rate of massive gauge quanta is exponentially suppressed. Gluon production is 
also suppressed unless a > Aqcd, that is, 

n/m > 1 GeV. (15) 

Thus we have for both photons and gluons with Q 2 < a 2 

e ~ (fi/em 2 ) 2 fit. (16) 

The production rate (per unit volume) of these particles is 

h i ^n M {t)e- l {dE/dt) il (17) 

where subscript i runs through 7 for photon and g for gluon. In particular for gluons we 
obtain 



le s (m 



(18) 



Let us assume that each gluon with energy e fragments to hadrons with a power-law energy 
spectrum KE~ q . If we take q = 1.5 the cascade limit will be somewhat weaker than in the 
case of a more realistic QCD fragmentation function. Let us proceed with this favourable 
case. Using the normalized fragmentation function 

N p (e, E)dE = ^^E~ 15 dE (19) 

and the rate of gluon production h g given by Eq. flTB"|) , one obtains for the diffuse proton flux 

1 , e s R p (E) (m 2 \ 4 ( Ey 15 dE 



Ip(E)dE =_ h? ^^) [t) _. (20 ) 

The cascade energy density is determined by the e-m radiation of monopoles and can be 
estimated as 



771^ 1 

u cas = 4vre 4 — -3. (21) 
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Requiring that u cas < u) b s , the proton flux at energy E ~ 1 • 10 11 GeV is bounded by 

E 3 I P (E) < 3.7 ■ 10 14 /i6 1 m~ W-W 2 , (22) 

where /i6 = /i/10 6 Ge^/ 2 . The symmetry breaking scale of strings is unlikely to be below the 
electroweak scale, /z > 10 4 GeV 2 , and the flux ( p2| ) is considerably lower than that observed. 



(iv) Necklaces. 



Necklaces are hybrid TD corresponding to the case N = 2 in Eq.(|lT|), i.e. to the case 
when each monopole is attached to two strings. This system resembles "ordinary" cosmic 
strings, except the strings look like necklaces with monopoles playing the role of beads. The 
evolution of necklaces depends strongly on the parameter 

r = m/ )id, (23) 

where d is the average separation between monopoles and antimonopoles along the strings. 
As it is argued in Ref. [j38|| , necklaces might evolve to configurations with r ^> 1, though 



numerical simulations are needed to confirm this conclusion. Monopoles and antimonopoles 
trapped in the necklaces inevitably annihilate in the end, producing first the heavy Higgs 
and gauge bosons (X-particles) and then hadrons. The rate of X-particle production is easy 
to estimate as |38] 

~ (24) 

Using Eq's(f|) and fl5|) one can calculate the fluxes of UHE protons and gammas taking 
into account the restriction due to cascade radiation 

uca. = \Ur^fj f Jvhji = l fy k (25) 
The separation between necklaces is given by |38| D ~ r -1 / 2 ^ for large r. Since r 2 fi 
is limited by cascade radiation, Eq.(^), one can obtain a lower limit on the separation D 
between necklaces as 

16 



D ~ f^-V > 10W10 6 GeK 2 ) 1 / 4 fcpc, (26) 
\ 4t w cas y 

where we used / w ~ 0.5. Another (weaker) constraint on the parameters of the model 
follows from the condition d > 5 S , where 5 S ~ l/(erj s ) is the string width and is the string 
symmetry breaking scale. This condition gives r max ~ r) m /r] s , where r\ m is the monopole 
symmetry breaking scale. For r\ m ~ 10 16 GeV and rj s ~ 10 3 GeV, one obtains r max ~ 10 13 , 
which corresponds to D min ~ 1 kpc. 

Thus, necklaces can give a realistic example of the case D « i? 7 , when the Universe is 
uniformly filled by the sources. The proton-induced EAS from necklaces strongly dominate 
over those induced by photons at all energies except E > 3 • 10 11 GeV (see Fig. 3), where 
photon-induced showers can comprise an appreciable fraction of the total rate. 

The spectra of protons and photons from necklaces are shown in Fig.3. The calcula- 
tions were performed with the help of Eq.([|) and Eq.([5]) using the SUSY-QCD fragmenta- 
tion functions. The dashed, dotted and solid lines correspond to the masses of X-particles 
10 14 GeV, 10 15 GeV and 10 16 GeV, respectively. The values of r 2 /i used to fit these curves 
to the data are 7.1 • 10 27 GeV 2 , 6.0 • 10 27 GeV 2 and 6.3 • 10 27 GeV 2 , respectively. They 
correspond to the cascade density u cas equal to 1.5 ■ 10~ 6 eV/cm 3 , 1.2 ■ 10~ 6 eV/cm 3 and 
1.3 ■ 10~ 6 eV/cm 3 , respectively, all less than the allowed cascade energy density for which 
we adopt the conservative value u cas = 2 ■ 10~ 6 eV/cm 3 . 

The absorption of gamma-radiation on the radio background is taken according with the 
upper and lower limits of fL9| , shown in Fig.2. The agreement with observational data at 
the highest energies is improved when the photon flux is added to the proton flux. 

Two following effects are expected to improve the agreement with the data at E > 
1 ■ 10 11 GeV. 

The fluctuations in energy losses increase the fluxes. This effect is significant when the 
fraction of energy lost in one collision becomes appreciable, i.e. at E > 1 • 10 11 GeV. This 
effect can be accounted for in Monte Carlo simulations. 
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The local enhancement of the density of sources, e.g. within Local Supercluster (LS), 
makes the spectrum at E > 3 • 10 10 GeV flatter (see calculations in ||56|| ). The accumulation 
of necklaces within LS is expected for large r. Indeed, for sufficiently large r, the typical 
velocity of necklaces v — cj ^Jr is less than the escape velocity from the Local Supercluster 
and thus necklaces are confined within LS. For r > 10 7 necklaces can also cluster on the 
galactic scale. 

For energy lower than 1 • 10 10 GeV we assume the presence of another component with a 
cutoff at E ~ 1 • 10 10 GeV . It can be generated, for example, by jets from AGN [[57], which 
naturally have a cutoff at this energy. 



(v) Monopolonium and SH relic particles. 



These two sources exhibit the same clustering property: they act as non-dissipative 
matter which clusters in the Universe in the same way as Cold Dark Matter. As a result the 
density of these particles in the Galactic halo is enhanced according to Eq.(|7]).0 The spectra 
of UHE protons and photons from decays of relic particles in the halo are calculated using 
Eq.(g) for the different distributions of X-particles in the halo given by Eq.(^|). The results 
are displayed in Fig. 4. The mass of X-particle is nix = 1 • 10 14 GeV. The solid, dotted and 
dashed curves correspond to the distribution (|9|) with (a, f3, 7) equal to (2,2,0), (2,3,0.2) 

2 Clustering of free monopoles in the galactic halo has been briefly discussed in Ref. |p6| . The 
authors argue that the monopole density hm is proportional to the dark matter density and that 
the X-particle production rate due to monopole-antimonopole annihilation is fix oc n 2 M . There 
are, however, some problems with this picture. The motion of free monopoles is strongly affected 
by the galactic magnetic field, and their density is not likely to follow that of CDM. Moreover, 
the probability for annihilation of free monopoles in the galactic halo is extremely small, and the 
resulting flux of cosmic rays is negligible for all reasonable values of the monopole density um- 
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and (1,3,1), respectively. The shape of the spectra is naturally the same for all three 
curves. Since we normalize the spectra by the value to fit the observational data at 
Ej ~ 10 10 GeV, the fluxes are the same in all three cases. The corresponding values of tiq 
are 6.3 • 10~ 42 cm~ 3 s _1 , 1.3 ■ 10 -41 cm~ 3 s~ 1 and 1.6 ■ 10~ 41 cm~ 3 s _1 , respectively. These 
values also normalize, through relation Eq.(|7|), the flux of extragalactic protons, given by 



Eq.(|]). The obtained results do not significantly differ from those in ||25|| . The photon flux 
is greater than the proton flux by a factor of 6. 

Another signature of this model is the anisotropy caused by asymmetric position of the 



Sun in the Galactic halo ||26|| . The anisotropy reveals itself most significantly as the large 
ratio of fluxes in the directions of the Galactic Center (GC) and the Galactic Anticenter 
(GA). We define the anisotropy A(9) as this ratio with the fluxes measured within a solid 
angle limited by angle 9 with the line connecting the Sun and the GC. This can be readily 
calculated using the distribution of the sources in the halo given by Eq.([|). 

The anisotropy A(9) is plotted in Fig. 5 as a function of 9 for all three density profiles 
given by Eq.(|9|) with r = 5 kpc (solid curves) and r = 10 kpc (dashed curves). The 
estimated anisotropy is relevant for energies E > 1 • 10 10 GeV, where the contribution of the 
isotropic lower-energy component is small. 

There is one more signature for this model, given by a direct flux of protons from the 
Virgo cluster, provided that they are weakly deflected by magnetic fields in the Local Su- 



percluster (see p5| ). Note that UHE photons from this source are absorbed, and cascade 



radiation does not propagate rectilinearly because of the cascade electron deflection in the 
magnetic field. 

(v) Vortons 

Vortons are charge and current carrying loops of superconducting string stabilized by 
their angular momentum Although classically stable, vortons decay by gradually losing 



charge carriers through quantum tunneling. Their lifetime, however, can be greater than 
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the present age of the universe, in which case the escaping X-particles will produce a flux of 
cosmic rays. The X-particle mass is set by the energy scale 7]x of string superconductivity. 

The number density of vortons formed in the early universe is rather uncertain. Accord- 
ing to the analysis in Ref. |58|, vortons are overproduced in models with rjx > 10 9 GeV, 
so all such models have to be ruled out. In that case, vortons cannot contribute to the 
flux of UHECR. However, an alternative analysis fl59f suggests that the excluded range is 
10 9 GeV < i]x < !0 l2 GeV, while for rjx 3> 10 12 GeV vorton formation is strongly suppressed. 
This allows a window for potentially interesting vorton densities withP] rjx ~ 10 12 — 10 13 GeV. 



Production of UHE particles by decaying vortons was studied in Ref. pOfl . As we already 
mentioned above, vortons cluster in the Galactic halo and the discussion in the preceeding 
subsection is therefore directly applicable to this well. 



V. DISCUSSION 

We studied observational constraints on various TD models, as well as possible signatures 
of TD as sources of the observed UHE radiation (E > 10 10 GeV). 

The most stringent constraint is due to electromagnetic cascades. It depends on the 
energy spectrum of particles from decays and on astrophysical quantities which determine 
the development of the cascade (most notably on the flux of intergalactic infra-red/optical 
radiation and on intergalactic magnetic field ). The SUSY-QCD spectrum makes the cascade 
constraint weaker, because this spectrum predicts less higher energy particles and more low 
energy particles as compared with ordinary QCD spectrum. In case of very large mx it 
means that for a given UHECR flux, less energy is transferred to the e-m cascade radiation. 
There are considerable uncertainties in the extragalactic flux of infra-red radiation and 



3 These numbers assume that strings are formed in a first-order phase transition and that r\x is 
comparable to the string symmetry breaking scale r\ s . For a second-order phase transition, the 
forbidden range widens and the allowed window moves towards higher energies p9|| . 
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extragalactic magnetic field. The conservative limit on the energy density of the cascade 
radiation imposed the latest EGRET data is uo cas m 2 • 1CT 6 eV/cm 3 . It could be 3 — 5 ■ 
10~ 6 eV/cm 3 with astrophysical uncertainties mentioned above. The further progress in the 
study of origin of EGRET extragalactic flux and in calculation of SUSY-QCD spectrum, 
in the pessimistic case, can exclude such TD as e.g. necklaces as the sources of observed 
UHECR. 

Another important constraint arises from the fact that at ultra-high energies, the proton 
attenuation length R P (E) and the photon absorption length R^(E) are both small compared 
to the Hubble radius. Models in which the typical distance between defects is D » R p are 
disfavored. In such models, the observed spectrum would have an exponential cutoff, unless 
a source is accidentally close to the observer. In the latter case the flux would be strongly 
anisotropic. 

Finally, in many cases TD give UHECR fluxes lower than the observed ones. We showed 
here that this is the case for monopole-string networks. Superconducting and ordinary 
cosmic strings probably belong to this category as well, although some loopholes still remain 
to be closed. 

With all these constraints taken into account, it appears that only necklaces, monopolo- 
nium and relic SH particles survive as potential UHE sources. 

The most important observational signature of TD as sources of UHE CR is the presence 
of photon- induced EAS. For all known mechanisms of UHE particle production the pions 
(and thus photons) dominate over nucleons. At energies lower than 1 • 10 12 GeV, protons 
have considerably larger attenuation length than photons and the observed proton flux 
can be dominant. Nevertheless, even in this case photons reach an observer from sources 
located inside the sphere of radius R 1 (E) (assuming that i? 7 > D). Unlike protons, photons 
propagate rectilinearly, indicating the direction to the sources. 

Necklaces with a large value of r = m/iid > 10 7 have a small separation D < i? 7 . 
They are characterized by a small fraction, R^/Rp, of photon-induced EAS at energies 
10 10 — 10 11 GeV. This fraction increases with energy and becomes considerable at the 
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highest energies. For smaller values of r ~ 10 4 — 10 6 , when the separation is larger than i? 7 
but still smaller than R p , most of UHE particles are expected to be protons (with a chance 
of incidental proximity of a source, seen as a direct gamma-ray source). Thus, in all cases 
necklaces are characterized by an excess of proton-induced showers. However, some fraction 
of photon- induced showers is always present, and it can be large at the highest energies. 

Monopolonium, decaying vortons and SH relic particles are characterized by an enhanced 
density in the Galactic halo. They give a photon-dominated flux without a GZK cutoff. Due 
to asymmetric position of the Earth in the Galaxy, this flux is anisotropic. The largest flux is 
expected from the direction of the Galactic Center, where the density of sources is the largest. 
Unfortunately, the Galactic Center is not seen by gigantic arrays , such as Akeno, Fly's Eye, 
Haverah Park, and Yakutsk array. However, these detectors can observe a minimum in the 
direction of the Galactic anticenter in particles with energies E > 1 • 10 10 GeV, as compared 
with the direction perpendicular to the Galactic Plane. 

A flux from the Virgo cluster might be another signature of this model. 

The search for photon induced showers is not an easy experimental task. It is known (see 
e.g. Ref. [51]) that in the UHE photon-induced showers the muon content is very similar 



to that in proton-induced showers. However, some difference in the muon content between 
these two cases is expected and may be used to distinguish between them observationally. 
A detailed analysis would be needed to determine this difference. 

The Landau-Pomeranchuk-Migdal (LPM) effect |62| and the absorption of photons in 



the geomagnetic field are two other important phenomena which affect the detection of 
UHE photons [|61|j63|| ; (see [21| for a recent discussion). The LPM effect reduces the cross- 



sections of electromagnetic interactions at very high energies. However, if the primary 
photon approaches the Earth in a direction characterized by a large perpendicular component 
of the geomagnetic field, the photon is likely to decay into electron and positron |6T| , |63| . Each 



of them emits a synchrotron photon, and as a result a bunch of photons strikes the Earth 
atmosphere. The LPM effect, which strongly depends on energy, is thus suppressed. If on the 
other hand a photon moves along the magnetic field, it does not decay, and LPM effect makes 
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shower development in the atmosphere very slow. At extremely high energies the maximum 
of the showers can be so close to the Earth surface that it becomes " unobservable" pl |. 



We suggest that for all energies above the GZK cutoff the showers be analyzed as can- 
didates for being induced by UHE photons, with the probability of photon splitting in the 
geomagnetic field determined form the observed direction of propagation, and with the LPM 
effect taken into account. The search for photon- induced showers can be especially effec- 
tive in the case of Fly's Eye detector which can measure the longitudinal development of 
EAS. The future Auger detector will have, probably, the highest potentiality to resolve this 
problem. 
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APPENDIX A: UHE PARTICLE PRODUCTION DUE TO MULTIPLE LOOP 

FRAGMENTATION 

We shall adopt the following simple model of loop fragmentation (It is somewhat similar 
to the model introduced in Ref. plfl ). Each loop fragments into (iVj + 1) daughter loops 
after it radiated away a fraction (1 — j\) of its energy. The Lorentz factor of the daughters 
in the center-of-mass frame of the parent loop is IY If the initial mass of the loop is M, 
then after n rounds of fragmentation, the number of daughters is 
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N n ~ N?, (Al) 
and their energy, rest mass and Lorentz factor are, respectively, 

E n ~ (fi/N^M, (A2) 

M n ~ (A/iV 1 r 1 ) n M, (A3) 

r n ~ r?. (A4) 

The fragmentation process stops at round n* when M n „ ~ r], where 77 is the symmetry 
breaking scale of strings, that is, when the size of the fragments becomes comparable to the 
string thickness, 

~ HwTnry (A5) 

Loops decaying at the present time t have masses M ~ 10 2 Gfi 2 t Q and 

M/77 ~ 10 55 ?7^ 6 , (A6) 

where rjiQ = i]/10 16 GeV . The typical values of j\, Ni and I\ are not known. Numerical 
simulations of loop fragmentation in Ref. |52| found N\ ~ 3 — 10 and Ti ~ 1.3. However, 



these simulations used loops of arbitrary shape, and the parameter values for loops in a 
realistic network may be quite different. It seems reasonable to assume that on average the 
loop has to lose at least 30% of its energy to initiate the next round of fragmentation, that 
is, fi < 0.7. As we shall see from Eq.([A9|), the energy output of loops in the form of cosmic 
rays is maximized for the largest possible values of N\ and j\. In the estimates below, we 
shall adopt Ni ~ 10, fi ~ 0.7 and I\ ~ 1.3. The values of Ni and fi appear somewhat large, 
and we shall keep in mind that we have made a rather optimistic choice of the parameters. 
With these values, 

ra* w 44 + In ?7i 6 . (A7) 
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Each fragmenting loop gives ~ N\ particles of energy ~ rjTi in the rest frame of the loop 
(~ 1 particle per intersection). The fraction F of the total energy of the initial loop that 
ends up in the form of UHE particles can be estimated as 

F ~T7lL N n?n. (A8) 
1V1 n=l 

The dominant contribution to the sum is given by the last term (that is, by the last round 
of fragmentation), and we can write 

^-^(Wirwr, (A9) 



where in the last step we have used the definition of n t . Using Eq. flA7j) and fi « 0.7, we 
have for the fraction of energy transferred to X-particles 

F ~ 2 • 10' V 6 °' 4 - (A10) 

Most of the particles are emitted at energies 

E x ~ r™*r/ ~ 9 • 10 20 7]% 7 GeV. (All) 

The X-particle injection rate is given by 

n x ~ F— ^-3 ~ 3 • W-^ri^cm-h- 1 . (A12) 

To maximize the diffuse proton flux, I P (E), we shall assume a power-law fragmentation func- 
tion KE~ P with p = 1.5 in the system where X-particle is at rest. After simple calculations 
using the Lorentz-transformation to the laboratory system, one obtains the diffuse flux as 

Ip(E) ^(l^lfK^(l.)-\ ( E), (A13) 



Anp Ex \Ex 

where R P (E) is the proton attenuation length. For E = 6.3 ■ 10 19 eV, using Eqs.( [AlI| , [A12j ) 



we obtain E 3 I P (E) w 3 ■ 10 18 ?^g 25 eV 2 mr 2 s~ l sr^ 1 to be compared with the observed value 
3 • 10 24 eV 2 m~ 2 s~ 1 sr~ 1 . Thus the calculated flux is too small by a factor of 10 6 . This 
discrepancy is difficult to resolve by stretching the range of the parameters fi and N±. For 
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example, with f\ ~ 0.9 (which appears unreasonably large), the proton flux is E 3 I P (E) « 
3 ■ 1 23 r7^g 5 eV/m 2 ssr, but the cascade energy density is too large due to increase of E x ■ 

u cas « {l/2)E x h x t w 7.4 ■ lO- 3 V li eV/cm 3 . 

We next consider particle production by cusp evaporation in the fragmenting loops. 
Assuming that cusps are periodically repeated and completely "evaporated" into particles, 
the energy rate of particle production (in the rest frame of the loop) is f3£| 



4 - M-r 1/3 . (am) 

V 



Compared to the gravitational radiation power, £ g ~ 10 G\i , this is 



S p /S g ~ lQ-\G^)-\—)-" A . (A15) 
The Lorentz factor at the cusp is 

r c ~ (^) 1 / 3 . (Ai6) 

The fraction of energy lost by a fragmenting loop in the form of particles is 

F = W-\GvT\l - h) £ )" 1/3 (A17) 

M V 



M ™ c 

= io- 2 ^)- 1 ^ - f^-y 1 ' 3 E(/i^ r i) n/3 » (Ai8) 

V Si 

where 

n c ~ 32 + 31n7? 16 (A19) 

is the value of n at which £ p ~ £ g . 



For our choice of parameters, f^N-^Ti > 1 and the dominant contribution to ( |A18|) is 
given by n ~ n c . Then 

F ~ (1 - /i)/r ~ 3 ■ 10-Ve 1 , (A20) 
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which is comparable to Eq. (|A10 ). The energy of the emitted particles is 



E x ~ riY^Y c ~ 4 ■ 10^6 GeV. (A21) 

This energy is too high: the observed flux cannot be obtained for any reasonable value of r\ 
without violating the cascade bound. 
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Figure Captions 



Fig. 1: Gamma-ray absorption length as a function of energy. The solid line and the 



dotted line reproduce the upper and lower case of ref. [19], while the dash-dotted line is 



taken from ref. [18] and drawn here for comparison. 



Fig. 2: The j/p ratio for different values of mx ( indicated in the figure) as a function 
of the gamma-ray energy. High (7-high) and low (7-low) photon fluxes correspond to two 
extreme cases of gamma-ray absorption from |19| . 



Fig. 3: Proton and gamma-ray fluxes from necklaces. High (7-high) and low (7-low) 
photon fluxes correspond to two extreme cases of gamma-ray absorption from ||19|| . The 
fluxes are given for mx = 1 ■ 10 14 GeV (dashed lines), mx = 1 ■ 10 15 GeV (dotted lines) and 
mx = 1 ■ 10 16 GeV (solid lines). The fluxes are normalized to the observed data. 

Fig. 4: Predicted fluxes from relic SH particles {mx = 1 • 10 14 GeV) or from monopolonia 
producing X-particles with the same masses: nucleons from the halo (curves labelled as "pro- 
tons"), gamma-rays from the halo (curves labelled "gammas") and extragalactic protons (as 
indicated). The solid, dotted and dashed curves correspond to (a, /3, 7) = (2, 2, 0), (2, 3, 0.2) 
and (1,3,1) respectively (see text). 

Fig. 5: Anisotropy A{9) as a function of the angle 9. A{9) is defined as a ratio of fluxes 
in the direction of the GC and AC. The fluxes are calculated within solid angles limited 
by angle 9 with the line connecting the Sun and GC. Anisotropy is given for three density 
profiles : (a,/3, 7) = (2,2,0) - "isothermal" curve, (2,3,0.2) - "best fit" curve, and (1,3,1) 
- for numerical simulation shape. The solid lines correspond to r = 5 kpc and the dashed 
lines to r = 10 kpc. 
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